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ABSTRACT 


It  is  shown  that  the  span  of  a  collection  of  multivariate  B-splines, 

S^  :=  span  M(*|x  °,*«*,x  k+m)  ,  contains  polynomials  of  total  degree  <  k 
jeJ 

if  the  index  set  {j  =  (jn,»»»,j,  ^  )} .  ,  has  a  certain  combinatorial 

u  K+m  jej 

structure,  in  particular  the  knots  {xV}  c  can  be  chosen  essentially 


arbitrarily.  Under  mild  assumptions  on  the  distribution  of  the  knots,  a  dual 
basis  jeJ  for  sk  *s  constructed  which  has  local  support,  i.e., 

supp  c  sUpp  ,  and  the  condition  number  of  the  B-spline  basis  is 
estimated.  This  leads  to  good  local  approximation  schemes. 
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Key  Words:  Multivariate,  B-splines,  spline  functions,  dual  basis, 
condition  number 
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SIGNIFICANCE  AND  EXPLANATION 


Spline  approximation  provides  good  approximation  methods  in  one 

variable.  This  is  to  a  great  deal  due  to  the  properties  of  the  B-spline 

0  n 

basis.  Recently,  multivariate  B-splines  M(*jx  ,**»,x  )  have  been  defined 
and  their  basic  recurrence  relations  derived.  However,  there  was  no 
satisfactory  way  to  choose  for  a  given  set  of  knots  {xV}  c  Rm  a  span  of  B- 
splines  that  contains  polynomials;  this  property  is  essential  for  good 
approximation  properties. 

In  this  report  we  give  a  sufficient  condition  in  terms  of  the  index  set 

j0  jjc+m 

{j  *  <j0»*** / jk  )}  .  j  for  S*  span  M(*|x  ,»*»,x  )  to  contain 

2  jeJ 

locally  polynomials  of  total  degree  <  k.  Since  the  knots  may  be  chosen 
almost  arbitrarily,  the  resulting  space  of  multivariate  splines  Sk  has  the 
local  flexibility  familar  from  the  univariate  theory.  We  show  that  under  mild 
assumptions  on  the  distribution  of  the  knots  the  L^-norm  of  a  spline  s  £ 
is  equivalent  to  the  i^-norm  of  its  normalized  B-spline  coefficients.  As  an 
application,  we  indicate  how  local  approximation  schemes  can  be  constructed. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MKC,  and  not  with  the  author  of  this  report. 


d 


JPPUfHf pp 


MULTIVARIATE  SPLINES 
Klaua  Hollig 

0.  Introduction 

'  The  purpose  of  this  report  is  to  introduce  a  class  of  multivariate  non 
tensor  product  spline  functions  and  study  its  basic  properties .^the  basis  for 
our  approach  are  recent  results  of  C.  A.  Micchelli  and  W.  Dahmen  [8-12,  21, 

22]  on  multivariate  B-splines  M(* |x°, • • • ,xn)  Introduced  by  C.  de  Boor  in 

,  ^0  Jn 

[3].  He  consider  a  space  of  spline  functions  S.  -  span  M(*|x  ,*»*,x  >  , 

J£J 

where  the  index  set  J  has  a  certain  combinatorial  structure  and  the  knots 
v  m 

x  e  R  may  be  chosen  almost  arbitrarily.  It  is  shown  that  contains 

polynomials  of  total  degree  k  «  n-m  which  is  the  basis  for  good  approxima¬ 
tion  properties.  Under  mild  assumptions  on  the  distribution  of  the  knots  we 
construct  a  dual  basis  A^,  j  e  J  ,  for  with  local  support,  i,e. 

supp  ^  c  supp  M ^  .  Using  extensions  of  the  functionals  A  ^  to  Lp  the 
condition  number  of  the  B-spllne  basis  is  estimated  which  turns  out  to  be 
essentially  independent  of  the  distribution  of  the  knots.  As  an  application 
we  obtain  the  basic  error  estimate  for  local  approximation  schemes. 

Since  the  spaces  have  the  local  flexibility  familiar  from  the 

univariate  theory  one  might  expect  applications  to  finite  element  methods, 
smoothing  of  data  and  surface  fitting. 


Sponsored  by  the  United  States  Army  under  Contract  No.  DAAU29-80-C-0041.  This 
material  is  based  upon  work  supported  by  the  National  Science  Foundation  under 
Grant  No.  MCS-7927062. 


1 .  B-Splines 


The  geometric  interpretation  of  the  univariate  B-spline  by  B.  B.  Curry  and  I.  J. 


Schoenberg  (7]  gave  rise  to  the  following  generalization. 


Definition  1.  C.  deBoor  [3].  Let  x°,**«,xn  be  any  points  in  tP,  n  >  m  ,  which  span  a 
proper  convex  set  and  choose  any  simplex  o< (x°,x°) , ••• , (xn,xn) )  with  vertices 
(xV,xV)  e  Rn  .  The  B-spline  M(4 |x°,*»* ,xn)iRra  +  R  corresponding  to  the  "knots"  xV  is 
defined  by 

n  vol  {x  £  o((x°,x°),***,(xn,xn))} 

.  _ .  ,  u  n  v  n-m 

(1)  M(x|x  .•••.X  ) - 5"=o - — - 

vol  0<(X  ,x  ),*»»,(X  ,X  )) 
n 


where  for  n*m  volQ(A) 


:■  ' 

.  i 


0  ,  if  A  -  d 


otherwise 


Denote  by  Afl  -  { (XQ ,• • • ,Xn) |  \  Xy  -  1,  >  0,  v  -  (>,•••, n}  the  standard  i 

n  v»0 

dimensional  simplex.  A  simple  calculation  shows  [22,  p.  3,4]  that  for  f  e  C(r“) 
_  .  t  ,  v,  „  _ r  ...  ..0 _ n, ....... 


(2)  nl  /  f(  l  -X  xV)dl  •••dX  -  /  M(x|x  ,• • • ,x")f (x)dx 

•  “  v  1  n 

v-o 

a  relation  used  by  C.  A.  Micchelli  to  define  the  multivariate  B-splineB  [21].  Equation  (2) 
defines  M  even  for  n  <  m  as  a  positive  measure  supported  in  the  convex  hull  of  the 
points  x°,»«*,xn.  From  (2)  the  Fourier  transform  of  M  can  be  computed  via  the  Hermite 
Genocchi  formula  for  divided  differences  [23,  p.  16] 


fln)(  I  xvV<V..dxn 


A  v«0 

n 


Applying  Paraeval's  identity  /fg  »  Jig  ,  hly)  ■  /  h<x)e  ixydx  ,  to  the  equation  (2)  we 


/  fi(x|x°,**» ,x")g(x)dx  -  /  M(x|xu,*»* ,xn)$(x)dx 
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A  repeated  application  of  ( 9 '  shows  that  the  distributional  derivatives 

DUM(«  |x°,***,xn),  |a|  ”  n-m*1  ,  being  linear  combinations  of  the  B-splines 
i,  1 

M( • | x  ) /  are  supported  in  convex  subsets  of  hyperplanes*  Hence  M  agrees  with 

a  polynomial  of  total  degree  <  n-m  on  every  subset  of  R®  not  cut  by  one  of  the  convex 
1 1  1 

sets  spanned  by  x  ,*»«,x  '  ,  0  <  i,  <  •••  <  i  <  n  [21,  22].  The  global  smoothness 

i  on 

y 

of  M  depends  on  the  geometric  configuration  of  the  points  x  .  If  every  subset  of  i 
points  of  {x°,*»*,xn}  forms  a  proper  convex  set  it  follows  again  from  (9)  that 

DaM  e  ,  |a|  -  n-Jt+1  .  Therefore, 

,  a  0  n.  ji™ [ 1 1  ^  id .  n— i .  it . 

(10)  M( * | X  ,***,X  )  e  (R  )  e  C  (R  )  . 

Perhaps  the  most  striking  identity  for  multivariate  B-splines  is  the  recurrence  relation, 
discovered  by  C.  A.  Micchelli  [21]  and  also  proved  by  W.  Dahmen  [B]  under  more  restrictive 
assumptions.  See  also  [16]  for  a  proof  via  Fourier  analysis. 


Theorem  2  [8,21,22].  Assume  that  the  points  x°,*#*,xn  e  Rm,  n  >  m  form  a  proper  convex 


set  and 

then  we  have 
(11) 


I  ^ 


v-0 


l  K  “  1  » 


V-0 


M(x |x  , •  •  •  ,x")  -  -2-  l  \  M(x|x°,*«»,xV-1,xV+  ,«»«,xn) 

n-  m  “  v 
v-0 

if  all  B-apllnea  occurring  in  this  equation  are  continuous  at  x  . 


Formula  (11)  expresses  M  recursively  as  a  convex  combination  of  positive  quantities 
and  thus  provides  a  stable  way  for  computing  multivariate  B-splines  [22].  So  far  we  have 
only  discussed  some  basic  results  on  multivariate  B-splines  and  we  refer  to  the  literature 
(8-12,  21,  22]  for  many  interesting  recent  developments. 
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2.  Some  topological  preliminaries 

We  need  some  basic  facts  about  triangulations  [14].  Since  there  seems  to  be  no 

convenient  reference  for  a  reader  with  little  background  in  topology  we  sketch  the  proofs. 

Let  (X. I )  denote  a  collection  of  nondegenerate  closed  slmplices  o, (X)  ■ 

i  i  1 

cr(x  °,».»,x  m),  i  -  (i  ,»<>*,i  )  e  I  c  a™*1  with  vertices  xv  e  X  c  am  .  (X.I)  is  called 

u  w 

a  triangulation  of  a  set  n  c  Rm  if  ft  «  u  a.  (X)  and  the  intersection  of  two  different 

i£I 

slmplices  is  either  empty  or  a  common  lower  dimensional  face.  Moreover  we  assume  that 
every  compact  subset  of  ft  is  Intersected  by  finitely  many  slmplices  only. 

Definition  2.  Denote  by  +  p^B  the  inscribed  balls  with  center  and  radius 

p^  for  the  slmplices  (^(X).  We  call  a  triangulation  regular,  if  there  exists  a  constant 
Y  such  that 

(12)  diam  ^(X)  <  y  .  lei. 

Lemma  1.  Let  (X,I)  be  a  triangulation  of  an  open  set  (l  c  if1  .  Then  for  any  finite  set 
of  vertices  {xv}v£j  c  x  there  is  an  e  >  0  such  that  (X,I)  with  xv  ■  xv,  u  i  J  is  also 
a  triangulation  of  ft  if  |xv  -  xv|  <  e  ,  v  £  CT  , 


Lemma  2.  Let  (X,I)  be  a  triangulation  of  a 
locally  finite  collection  of  slmplices  (X,I) 
also  a  triangulation  of  ft  . 


simply  connected  set  ft  and  assume  that  a 

satisfies  ft  «  u  c.(X).  Then  (X,I)  is 
iel 


Proof.  The  mapping  f  t  ft  ♦  0  defined  by 


mi  v»  '•v 
fix  )  »  x  , 


f,  i  i8  «ffine  linear, 

I Oj i") 


is  a  local  homeomorphism  which  is  by  assumption  surjective.  Since  ft  is  simply  connected 
the  monodromy  theorem  implies  that  f  is  bijective  and  hence  (X,I)  is  a  triangulation  of 
ft  .  Roughly  speaking,  the  argument  is  as  follows.  Suppose  there  exist  £  3, 

•  *nch  that  y  ■  f(zQ)  ”  ftz^)  and  let  Cn  be  the  image  of  a  curve  joining  zQ 
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and  Zj.  Danote  by  C(t,X):  [ 0 , 1  ] 2  ♦  (3  a  continuous  deformation  of  Cg  «  C(»,0)  into 
tha  point  y  ,  i.e.  y  s  C(», 1).  Since  we  have  assumed  that  only  finitely  many  simplices 
intersect  a  given  compact  set  every  point  x  e  Cl  has  a  neighborhood  CJX  such  that 
f  'Ojj  consists  of  finitely  many  disjoint  sets  and  f(v  is  a 

_  1  X/  V 

homeomorphism.  Hence  the  inverse  of  f  along  each  curve  C^»  C( • , A) ,  f  '(C( •,!)),  with 
initial  value  zQ  -  f_1(C(0,X))  is  well  defined.  We  shall  show  that  -  f_1(C(1,X)) 

for  all  X  c  [0,1]  contradicting  f  \C(«,  1))  =  z^.  Por  a  sufficiently  fine  subdivision 
of  [0,1]  into  squares  Ry  each  of  the  sets  C( R^j  is  contained  in  some  neighborhood 
Ux  and  the  following  picture  shows  that  inverting  f  along  adjacent  curves  leads  to  the 

same  value  , 

A 

X 

f-1(C(i.A))  -  f“\C(t<1)>)  -  f_1(C(l,X’)). 


The  combinatorial  product  of  the  sets  i-  {i. ,•••,!  }  and  {0, •  ..,k}  is  defined  as 

o  rn 

the  collection  of  all  ordered  subsets  of  i  X  x  {0, • • .,k}  with  n+1  s-  m+k+1 

0  in 

distinct  elements,  i.e. 


(13) 


iAJc  s-  {( j,r)  -  ( 


r0 '  ***'rn 


)  e  •••,  im>  x  {0 ,  •  •  • ,  K>  ] 


n+1 


(3v+1  “  K  4nd  rv+1  ■  rv+1)  °r  <3v+1  *  V1  an<i  rv+1  “  V  for  v  ■  0' ••••")• 

Note  that  (j,r)  e  iik  is  already  uniquely  determined  by  j  -  (j  ,  ...,j  )  so  that  we  may 

0  n 

drop  r  when  referring  to  elements  of  ifik  . 


!■— a  3  .  Denote  by  ^  “  o(*0>  <  with  eQ  -  (0,  •••,()),  e1  «  (1,0,  •••,0),  •••,  the 

standard  k-dlraensional  simplex  and  let  (X,I)  be  a  triangulation  of  ,1  c  r"  .  Then  the 
simplices 
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3»  Splines  as  linear  combinations  of  B-splines 
Denote  by  the  polynomials  of  total  degree  <  k  ,  i.e. 

m  m  a 

P  *  {f|D  f  =*  0,  |al  =  k+1 } ,  where  a  *  ),  lal  *  J  a  ,  Da  *  II  (~“)  V  • 

K  i  m  v  -  3y 

v*  1  v 

We  call  a  space  of  splines  of  total  degree  k  defined  on  a  set  ft  c  Rm  if 

(A)  is  the  linear  span  of  a  collection  of  B-splines  and 

<B>  Pk(0)  c  Sk|n  . 

Moreover  we  shall  assume  that 

(C)  every  compact  subset  of  £5  is  intersected  by  the  supports  of  finitely  many 

B-splines  only. 

In  the  univariate  c*,e  we  may  associate  with  any  increasing  sequence  of  knots  {t^}, 

t\i  <  fc\rt-k+1  '  a  3pace  of  splines 

(15)  S.  =  span  M(  •  1 1  .••*,t  , , 

k  u  v+k+1 

v 

In  several  variables  any  triangulation  ((X,X),Z)  of  a  product  ft  x  ft  c  nm  x  rk, 

vol  ft  <  »  ,  corresponds  by  Definition  1  to  a  family  of  B-splines  which  forms  a  partition 
of  unity  over  ft  [3],  i.e. 


~  -1  -A  i  i  i  i 

(16)  1  z  (vol  ft)  £  vol  o((x  ,x  (x  n,x  n))  M( x | x  n)  , 

i*! 


x  e  ft,  n  *  m+k. 


and  it  was  shown  in  [91  that  for  ft  x  ft  =  [0,1]n  c  span  M(»|x  ,  •■•,x  n)  ,  i.e. 

iel 

condition  (B)  is  satisfied. 

To  require  the  existence  of  an  associated  triangulation  is,  however,  too  restrictive 
and  not  quite  satisfactory,  since  the  B-splines  do  only  depend  on  the  knots  xV  t  X  c  Rm  . 
As  in  the  univariate  case,  one  would  like  to  associate  with  any  distribution  of  knots  a 
corresponding  space  of  splines.  This  can  be  accomplished  by  a  multivariate  analogon  of  the 
process  of  "pulling  apart"  knots  which  we  shall  now  describe. 
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1 


Definition  3.  Let  (X,I), 


{xV}. 


,  be  a  triangulation  of  tP  and 


(xV)  xv  £  *“ 

1  lTv£X,u«0,***,k  '  »  ' 


a  collection  of  points,  not  necessarily  related  to  (X,I).  He  define  the  space  of 
splines  of  total  degree  k  with  respect  to  the  index  set  I  and  the  knots  X  by  ( c.f . 


(13),  (14)) 
(17) 


Sk,X,I 


span  M. (• |x) 
j£lik  3 


where  M^(«|X) 
satisfied. 


3n 


M(  •  |x 


■••• ,xr  ), 
n 


n  -  art-k 


In  addition  we  assume  that  condition  (C)  is 


The  rest  of  this  chapter  is  devoted  to  the  justification  of  this  Definition,  i.e.  to 
prove  (B).  To  simplify  notation  we  shall  drop  most  of  the  subscripts  in  the  sequel,  e.g. 
Sk  i“  “j  «"  ( • |  Z) ,  ai  j-  oA(x),  etc. 

Lemma  4.  if  we  make  for  the  knots  X  the  particular  choice  x^  -  xV,  v  £  *,  u  “  0,  •  • • ,k, 
then  j.  coincides  with  the  space  of  piecewise  polynomials  of  total  degree  k  with 
respect  to  the  triangulation  (X,I),  i.e. 

(18)  Sk,X,I  "  Pk,X,I  **  {f|f|0i(X)  6  Pk  '  1  £  I}  * 


V  V 

Proof.  For  the  special  choice  x  ■  x  we  have 

-  ^  p 


In 


k,X,l|a, 


“  span  M( • | x  ,»»»,x  ) 

j£i4k 


and  (x  .♦•*,x  "}  “  {x  m]  .  We  could  new  use  the  formula  for  B-splines  with 

multiple  knots  (22,  p.  16)  to  complete  the  proof.  But  we  may  also  apply  (9)  to  conclude 
that  the  supports  of  the  distributional  derivatives  D°M,  I u |  ■  k+1,  are  contained  in  the 


boundary  of  the  simplex  a  .  This  implies  span  M  c  p  (o,).  The  reverse  inclusion  is  a 

j£  iAk  3  *  ‘ 

special  case  of  Theorem  3  below.  Moreover,  since  #iAk  •  (  ^  )  ”  dim  P^,  the  B-splines 
My  j  £  iAk  ,  form  a  basis  for  P^to^. 
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Remark .  Lemma  4  shows  that  we  may  view  the  spaces  as  perturbations  of  piecewise 

v 

polynomials  obtained  by  pulling  apart  the  multiple  knots  x  e  X 

u  u  „ 

x  ♦  xy  ,  u  -  0,**»,k  . 

m  V  V 

One  may  e.g.  choose  a0,**«,ak  £  R  and  define  x^  ■  x  +  ay  ,  u  »  0,«*»,k.  This  very 
special  perturbation  corresponds  to  an  affine  transformation  of  the  triangulation  (14)  and 
is  only  an  appropriate  choice  if  the  triangulation  (X,I)  is  quasiuniform  and  one  is 
interested  in  a  uniform  distribution  of  the  knots  X  . 

We  have  defined  the  spline  spaces  on  R01  to  avoid  a  separate  discussion  of 

various  types  of  "boundary"  conditions  in  the  sequel.  This  is  no  loss  of  generality  since 
we  include  all  kinds  of  finite  dimensional  spline  spaces  by  simply  restricting  the  domain 
of  definition.  The  global  smoothness  of  depends  on  the  smoothness  of  the  particular 

B-splines  (c.f.  (10)).  But,  although  the  set  of  knots  yielding  the  highest  possible 
smoothness  c  C*-1  is  dense,  there  is  no  canonical  way  to  choose  such  knots  as  in  the 
univariate  case. 


Example  1 .  Let  {ty} ,  ty  <  tv+k+1»  be  an  increasing  sequence  of  knots.  To  obtain  the 

canonical  B-spline  basis  (15)  we  set  (c.f.  Definition  3)  X  ■»  {t,,.,.  }  _  , 

(k+1  )v  ve* 

I  -  {<v,v+1)}yeZ  ,  xjj  «  ‘(k+l  )v-u  '  v  E  *'  u  "  Which  yields 


(v,v+1  )Ak 


Kv«  •  «v(v+1 ) 

0  •  •  •  p  p 


(v+1h  „ 

k  >'  p 


0 


.k}  , 


span 

j£(v,v+1  )Ak 


M. 

3 


span 

*0, • • • , k 


M(*  |t 


(k-H)v-p 


.  •  •  •  •  t  )  « 

lk+1 ) (v+1 )-p ' 


The  following  example  shows  that  Definition  3  is  slightly  more  general  than  the  usual 
univariate  definition,  since  we  do  not  require  that  a  B-spline  does  only  correspond  to 
adjacent  knots. 


Example  2.  Let  X  -  Z,  i.e. 


v, 


{ (v,v+1 )} 


veZ 


and  consider  the  knots 


2, 


1,  x. 


1 


2.5,  x* 


v 

0.5  and  x 

V 


v  for  v  *  1,2 


Restricting  the  domain  to  the 


interval  [0,3]  we  have 


S,  ~  rim  ,,  “  span{M(*  lo, 0,2.5),  M(*  |0,2,2.5),  M(  •  |0.5,2, 2.5) ,  M(  •  |  0.5, 1,2 ) , 

1  t  A  I  l  U  t 

M(* | 0.5, 1,3) ,  M(« | 1,3,3)} 

which  is  a  space  of  continuous  piecewise  linear  splines  with  an  unusual  basis. 


Theorem  3.  The  spaces  x  1  contaj-n  polynomials  of  total  decree  k  .  In  particular 
have  for  all  x,y  e  R® 

09)  ( 1+xy)*  -  l  C4(y|X)  M„(x|X) 

where  (c.f.  (13),  (14),  (17)) 


jelAk  J  3 


C  (y| X)  i«  e  det 
j  3  nl 


and  e^  €  {-1,1}  is  chosen  so  that 


C.(X)  s-  e.  77  det 
3  j  nl 


30 

d  ♦  *  ®y,e 
0  0 


(1  ♦  xS.r 
n  n 


e  ••• 


1  •••  1 


>  0 


(Recall  that  xW  are  the  vertices  of  the  trlangulatlon  corresponding  to  the  index  set  I)  . 

Por  the  standard  univariate  B-spline  basis  formula  (19)  is  known  as  Marsden's  identity 

[19].  The  basic  idea  of  its  multivariate  generalization  is  due  to  H.  Dahmen  [9].  For  two 

variables  (19)  was  proved  by  T.  N.  T.  Goodman  and  S,  L.  Lee  [15]  under  the  restrictive 

2 

assumption  that  the  B-splines  correspond  to  a  triangulation  of  R  x  in  both  cases 

there  is  a  more  explicit  formula  for  the  coefficients  C^(y|X)  (c.f.  Corollaries  1  and 


2).  Besides  the  generalization  to  an  arbitrary  number  of  variables  we  show  that  identity 


9 


m 


(19)  is  valid  for  any  collection  of  knots  as  long  as  we  keep  the  costal natorial  structure 
determined  by  the  index  set  1  . 

Remarks .  (i)  Identity  (19)  holds  for  any  finite  collection  of  B-splines  M ^ ,  j  t  j  c 

IAk,  if  we  restrict  x  to  the  domain  Rm  \  U  supp  M  . 

j/J  3 

(ii)  For  y  ■  0  ,  (19)  is  a  generalization  of  (16)  since  the  coefficients  i»  C(0|X) 
need  not  to  be  positive.  The  coefficients  corresponding  to  the  B-splines  in  Example  2, 
e.g.,  are  1.25,  1,  -0.8,  -0.4,  1.25,  1.  One  can  show,  using  Lemma  2,  that  the 
coefficients  C.,  j  £I&k,  are  positive  if  and  only  if  the  simplices 

3  3  3 

(20)  o,(X)  s-  0((X  °,e  ) ,  *»*,(x  n,e  )),  )  £  Uk  , 

3  0  0  n  n 

form  a  triangulation  of  R  x  A  . 

(iii)  Comparing  the  coefficients  of  ya  in  identity  (19)  we  obtain  an  explicit 
representation  for  the  monomials  xa,  |a|  <  k  , 

(21)  x°  -  l  c“(X)  M  (x) 

3  fidk  3  3 

where  C?(X)  s-  (al f1  <  ,*,  )"3  D°C  (y|X|.  .  ,  a!  -  H  at  . 

3  lot  3  ly-o  ^  v 


Proof  of  Theorem  3.  Note,  that  by  Definition  1  we  may  associate  with  the  B-splines 
the  simplices  a ^  (c.f.  (17),  (20)).  First  we  assume  that  the  collection  of  simplices 

Oy  3  «  lAk,  is  a  small  perturbation  of  the  triangulation  (14),  l.e.  x  j  is  close  to 
pk,X,I  <c*f*  (18)).  More  precisely,  for  a  finite  set  J  c  z  and  e  >  0  ,  to  be  chosen 
later,  we  assume 


<•) 


(*> 


I -  xv|  <  e,  u  -  0,..«,k,  v  e  J 


(**)  xj|  ■  xv  ,  u  ■  (),•••, k,  v  f  J  . 

By  Lemma  1,  ,  3  lAk  ,  is  a  triangulation  for  small  e  (Since  Rm  x  c  Rn  is  not 

an  open  set.  Lemma  1  cannot  be  applied  directly.  But,  we  may  consider  suitable  extensions 
of  the  triangulations  to  I?1).  Let  K  be  an  arbitrary  bounded  subset  of  Hf”.  Obviously  it 
is  sufficient  to  prove  (19)  under  the  assumption  x  r.  K  and  y  small,  vfe  define  the 
mapping 
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_*  _k  _■  _k 

T  i  K  x  X  ♦  *  x  ft 

T(x,x)  »  (x,(1  ♦  xy)x)  . 

Consider  the  finitely  many  B- splines  ,  j  <  j  ,  tha  supports  of  which  intersect  K 

(c.f.  Cha  assumption  in  Definition  3),  and  let  (xV)  „  .  „  ,  ba  tha  corraa ponding 

U  V€J 

aat  of  knot a.  ror  aaali  enough  y  tha  simpliees 

T  3o  ^n  - 

o*  i“  o (T(x  #e  ),*••, T(x, a  )),  i  t  J  , 

3  0  0  n  n 

font  alao  a  triangulation •  Moreover,  tinea  T  nape  tha  hyparplanaa  bounding  tP  *  into 

hyparpianaa  and  laavaa  pro j act Iona  on  rf1  fixed  wa  have 

t  4.  o’]  n  X  x  »k  -  T(K  x  A  )  . 

JeJ  3 

By  Definition  1  and  (16)  it  followa  that  for  «Mll  y  and  x  e  X 

vollx  £  W'Kx.x)  £  T(X  X  Ak»  -  (kl)“\l  +  xy)k  - 

VOl{x  f  Mk|(x,x)  £  0*}  -  vol  0*  M  . 
j  £J  3  j£j  3  3 

ror  Ix^  -  x^|  ,  v  £  J,  aaali  wa  have  agn  C^(y|X)  ■  agn  C^(X)  “  1  which  lapliaa 

T 

vol  c^  a  (kl)  C^(y|X)  and  coaiplataa  tha  proof  under  tha  aaauaptlona  (*),  (••). 

Ha  now  drop  tha  aaaiuxptlon  (*).  Let  11^,  j  £  J,  ba  tha  finitely  aany  B-apllnaa  having 
at  laaat  one  of  tha  knota  {x,V}  .  „  ,  and  consider  tha  aquation 

M  M™v  g  9  •  *  fK  t  V  €  J 

I  <1+xy)Vx)dx  -  /  (  J  C  <y|X)M.(x))6<x)dx  - 

tr  f  j£lAk  3  3 


J  I,  c  (y|X)M  (x)6(x)dx 

tr  j£j  3  3 


where  )  ,  )  c  i)  la  tha  Fourier  tranafora  of  a  teat  function  with  cxxapact  support.  By  (4) 

and  Paraaval'a  identity,  the  right  hand  aide  ia  equal  to 

nl  [  C  (y|X)  /  (-1*  °x,*”,-ix  "xl  •t  *(x)dx  . 
j£J  3  *  0  n 


The  coefflclenta  C^(y|X)  and  tha  divided  dlffarance  of  tha  exponential  function  are  also 

defined  for  complex  values  x^  £  C**,  and  it  is  easily  seen  that  tha  right  hand 

r  v«  m  U*  ( k+ 1 ) 

side  of  (a)  ia  an  entire  function  of  tha  arguments  lx  }  ,  £  (C  ) 

Tha  left  hand  aide  of  (+)  la  constant  with  respect  to  these  arguments.  Therefore,  since  by 
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the  first  part  of  the  proof  th«  identity  is  valid  in  a  raal  neighborhood  of  the  knots 

v 

x  ,  v  t  J  ,  it  follows  globally  by  the  identity  theorem  for  power  series  (17).  To 
obtain  (19)  frosi  the  identity  (+)  we  note  that  {+|+  t  V)  is  dense  in  V  with  respect  to 
the  norsi  Ifl  »•  sup  (1  +  |x|)4|(i(x)|  for  all  id. 

The  general  case  follows  now  easily,  consldsr  an  arbitrary  perturbation  of  the 
initial  knots 

V 

x  *  x^  ,  p"0 , •  •  •  ,k 

such  that  any  bounded  subset  It  c  ^  is  lnt.ersactsd  by  the  supports  of  finitely  many  b- 
splines  only.  Restricting  x  to  K  ,  only  finitely  many  knots  are  relevant  for  the 
Identity  (19)  and  the  previous  arguments  apply. 

we  now  examine  in  more  detail  the  coefficients  C^(y|S).  Indicating  torus  of  the  fora 
(1  ♦  *;'yi  by  "  ' s",  the  determinant  in  their  definition  has  the  following  structure 


1  ...  1 

where  the  poaltlon  of  the  nonsero  entries  in  the  first  k  rows  depends  on  the  index  set 

r"<r„,..*,r).  If  p  +1  ,  v»  ),•••, k,  denotes  the  number  of  's  in  the  v's  row 
On  v 

we  have  0q+“**(\  “  "  *  Por  suall  ■  this  loads  to  considerable  siupllf lcatlon. 

Corollary  1  119].  for  any  increasing  sequence  of  knots  {t  )  c  R  ,  t  <  t^^^j  ,  we 

have 


(23) 


(y-x) 


(k*  1 ) 


-1 


(t 


u+k+1 


-  t  ) 


k 

1  (y- 
u “) 


t  )  M(x|t  , 


•  •  •  1 1  .  .  ) 


Corollary  2  [15].  Let  be  defined  on  **.  Then  we  have  for  all  x,y  t  l2 

(24)  (1+xy)k-  l  C  a  (l+s^y)  M,(x) 

v  v  3  v-1  33 

where  x^  «•  s^(X)  can  toe  determined  by  the  knots  of  the  B-spline  Mj  and  the  index 
aet  j  •  For  details  we  refer  to  (IS]. 


Both  Corollaries  can  be  obtained  from  Theorem  3  by  direct  computation  of  the 
determinant  (22).  Mote  that  in  (23)  we  have  made  the  substitution  y  ♦  -1/y  • 
Unfortunately,  an  identity  of  the  form 

k  w 

C  (y|X)  -  C  a  (1  +  Xj (X)y) 

3  3  y.,  3 

ie  no  longer  valid  in  more  than  two  dimensiona. 


o  (ttyj-Yj)  (t-*y2-*y3)  o  o 

0  0  o  (i-*y1-yJ>  d+  y2*y3> 

o  i  o  t  o 

0  0  10  1 
i-i  i-i  i 

ii  ii  i 

cannot  be  written  in  the  form 

C<1  ♦  s,y) ( 1  ♦  s2y) 

for  same  vectors  s2  e  B3. 
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The  standard  univariate  B-spline  basis  (IS)  has  the  striking  property  (c.f.  C.  de  Boor 
[1,3])  that  the  Lp-norm  of  a  spline  is  equivalent  to  the  norm  of  its  coefficients  with 

respect  to  a  normalised  B-spline  basis  independently  of  the  sequence  of  knots.  Under  mild 
assumptions  on  the  distribution  of  knots  we  shall  establish  this  result  for  the 
multivariate  spline  spaces  Sv  _  T. 

Definition  4.  He  call  X  •  {x^}  a  regular  perturbation  of  a  triangulation  (X,I)  if  there 
exists  a  constant  x  and  balls  t-  +  hj.B  such  that  for  every  i  t  I  (c.f.  (13), 
Definition  3) 

(25)  0  c  n  supp  M 

j«iAk  3 

(26)  «i  n  supp  -  4  ,  j  e  l'Ak  ,  i>  *  1 

(27)  diam  supp  Mj  *  *  '  j  e  iAk 

Note,  that  the  conditions  (25)- (27)  are  local,  i.e.  we  do  not  assume  any  relationship 
between  hj,  hj,,  i  *  i’  . 

The  following  example  shows  that  assumptions  (25)-(27)  cannot  be  essentially  weakened. 
Example  4  Consider  the  following  perturbation  of  the  standard  triangulation  (14)  of 
X  *  [0,1]  restricted  to  [0,3]  *  [0,1] 


0  12  3 


Among  the  corresponding  B-spllnes,  M(*|l,2,2.5)  occurs  twice,  a  possibility  not  excluded 
by  Definition  3.  Although  the  B-splines  became  linearly  independent  if  we  move  the  vertex 
(2,1)  slightly  totheleft,  the  basis  is  not  well  conditioned. 


-T  -  "  n  " m  —  ’  - 


lemma  5»  With  the  notation  of  Definitions  2-4  unm  that  xi  can  be  chosen  as  the  center 
of  the  Inscribed  ball  of  the  ■  implex  o  ( x ) ,  i.e.  a  ,  and 

(28)  |x^  -  xV|  <  hv 

where  hy  «-  Min{pi  -  hi|xW  is  a  vertex  of  o^X)).  Then  X  is  a  regular  perturbation 
of  (X,I). 

This  shows  that  conditions  (25)- (2?)  allow  perturbations  of  the  same  magnitude  as  the 
perturbation  of  Example  4. 

Proof.  The  support  of  each  B-spline  M , ,  j  e iAk,  is  the  union  of  slmpllces 

““ 11  i  i  3 

o^  i»  o(x  °, •••,»!  *)  where  |x  V  -  x  W|  <  -  h^  .  By  T emma  2  the  collection  of 

siaplices  is  a  triangulation.  To  prove  (25),  (26)  it  is  sufficient  to  show  that 

for  any  such  triangulation  8^  <=  o^.  We  define  the  siaplices 

-_i  1  ^1  i 

«.  ■  o(tx  ♦  ( 1— t)x  ,***,tx  "  ♦  (1-t)x  "),  t  e  (0,1) 

*  »V 

and  observe  that  none  of  the  (k-1 {-dimensional  faces  of  a.  ,  t  e  [0,1],  intersects 

i  ,t 

0  0 

Sj  .  This  follows  because  any  such  face  can  be  separated  frca  8^  by  the  hyperplane 
parallel  to  the  corresponding  face  of  O^IX)  touching  8^  .  Concerning  (27)  we  note  that 
under  the  assumptions  of  the  Lena*  <  <  YP1+Pi~hi  . 

Theorem  ♦.  Let  X  be  a  regular  perturbation  of  a  trianguletlon  (X,l).  With  the  notation 
of  Theorem  3  and  Definition  4  we  define  functionals  1  i  C*t(rf')  ♦  ■  by 

(29)  X  (X,T)f  •  l  (ol)_2(  ,k.  )*’  d“c  (0|X  -  T)  D°f(T),  j  £  lit  , 

3  |u|<k  |a|  3 

where  T  £  8^  end  X  -  t  i -  { x|| — x)  is  a  translation  of  the  point  set  X.  Then  for  any  set 
of  points  t ^  £  8^,  j  £  14k  ,  the  functionals  X^  i»  X^(X.t^),  j  £  IAk  ,  are  dual  to 
the  B-sollne  basis  for  8^  !••• 

(30)  XJ  WJt  -  «jJt  ,  J.j’  £  IAk. 
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He  note  the  following  immediate  consequences  of  this  result  which  for  univariate 


spline  functions  is  due  to  C.  de  Boor  [2] . 

(i)  The  B-splines  M.  ,  j  e  IAk  ,  are  linearly  independent  over  the  set  ft  -  U  n  . 

3  id  1 

(ii)  For  any  polynomial  of  total  degree  <  k,  p  e  (X, tip  is  constant  as  a  function 

_  • 
of  T  £  R  . 


To  see  this  we  observe  that  the  polynomial  g(  x)  :»  X^ (X,  x)p  is  constant  for 
t  £  ft^  j  £  ifik.  This  follows  from  (30)  since  by  (26)  and  Theorem  3 


(31) 


sk.x.i|ni  '  ^  "ila^  ^ 


Concerning  (i)  we  note,  that  in  contrast  to  the  univariate  case  in  more  than  one  variable 
the  B-spline  basis  of  is  not  linearly  independent  over  arbitrary  open  subsets  of  tP 


Example  5.  Consider  the  following  set  of  knots  in  R2 


♦ 


which  may  be  extended  to  a  regular  perturbation  of  a  triangulation  of  R2  and  the 

,0012  ,0112 

corresponding  piecewise  linear  B-splines  N( •|Xq,x^,x1,x^ ),  M( • |x0,xQ ,x1 ,x1 ) , 

M(  •  |  xjj  ,x3  ,Xg  ,x^ ) .  Since  also  supports  of  other  B-splines  overlap  the  set  ft,  ft  contains 
a  set  of  linear  dependence  ft  ,  i.e.  the  B-splines  {H^lsupp  n  n  *  4}  are  linearly 
dependent  over  ft  . 


Proof  of  Theorem  4.  Since  supp  =  ft^  ,  j  £  iik  ,  and  S^j^  ”  Pk|ft  we  'iave  to  S*K>W 


M  l  »vMV>  “ 

3  j'eidk  3  3  3 


for  all 


polynomials  £  a  c  P.  •  A  change  of  variables  in  identity  (19) 

J’eiAk  3  3  ,ni  ,Sii 


*  *» 


yields 
(19') 

Wo  choose  T  ”  Tj  '  aPP1l'  the  functional  to  the  left  hand  side  and  obtain 


(1  +  <x-T)(y-T))k  -  l  c  ,  (y-T  |X-T  )M  ,  (x) ,  xefj  . 

j’elAk  3 

>  apply  the  functional  A^  to  the  left  hand  side  and  ot 

M  l  <,»,  )(x-T  )B(y-T  )6)  -  l  (al)_1D°C  (0|*-T  ){y-T  )“ 
3  |<5|<k  101  3  3  |a|<k  3  3  3 


Since  the  set  of  polynomials  {(1+(x-t  )(y-r))  }  spans  P.  ,  the  proof  is  complete  , 

J  j  K 

J  yeR 

if  we  show 

I  (a i )-1DaC  (0|X-T . ) (y-T . )a  -  C,(y-T. |X-t.),  y  £  r“  . 

|a|<k  3  33  3  3  3 

But  this  follows  simply  by  differentiating  with  respect  to  y  at  y  ■  . 

Formula  (29)  is  only  one  particular  representation  of  the  dual  basis  for  Mj,  j  e  iAk. 

We  now  construct  bounded  extensions  A_.  :»  A^(X)  of  the  functionals  A^  to 

Lp  ,  i.e.  AjtLp  -  R 

(32)  A^s  «  A^s  ,  s  £  . 

If  we  require  supp  A^  c  (3^,  j  £  iAk,  the  system  (32)  reduces  to 


(33) 


V  “  V  '  *  £  Pk  • 


Hence  we  may  determine  A^  from  the  identities  (19),  (21).  After  a  change  of  variables  we 
have 
(21*) 


I  (ol  )_1  ( |\  )_1  D°C,  (0|X-i)M  (x) 
jelAk  |a|  3  3 

and  clearly,  (33)  is  equivalent  to  the  system 


(33’) 


Aj (•  -  T)a  -  (0U)-1(|£(rV,Cj(0|X-T),  |a|  <  k  . 


Let  Lq  be  any  set  of  functions  dual  to  the  monomials  on  the  unit  ball  B  ,  i.e. 


(34) 


/  Lo<x)  xBdx  -  «ag  ,  | a | ,  |Bl  <  k  . 


We  transform  the  functions  h  to  the  balls  (J.  (c.f.  Definition  4)  and  normalize  them 

01  i 


in  Lpi 


(35) 


♦a(x)  -  h”nv/p'  L<j(h",(x-x1)) 


'Vp-.Oj  * 


i.e.  *  represents  a  bounded  functional  on  L  ((!.  )  .  we  now  set 
a  P* 


(36) 


[  *  ♦. 
■  a  a 


|a|<k 
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<*-  » 


and  determine  the  coefficients  from  03' )•  Choosing  t  •  x  ,  j  e  iAk  ,  we  obtain 


(37 ) 

which  implies 

(38) 


A  .  (•  -  x. ) 


I  a  h  "  p  /  L  (x)(h  x)“dx  -  a  h“ 
|Bl<k  P  1  b  @  1  “  1 


m/p*| a  I 


aa  “  <a,)"1(k|)"1hI^P"|0t  ^(OIX  -  x.) 


To  simplify  statement  and  proof  of  the  theorem  below  we  introduce  the  following 
notation.  With  a  space  of  splines  St  „  T  we  associate  a  partition  U  «  Rm  such  that 

for  each  j  c  IAk  there  exists  a  set  of  indices  Jj  such  that 

(39)  supp  M  -  u  U  . 

3  V£j 

Moreover  we  define  3 

(40)  Jy  i-  {J  £  IAk | supp  n  uv  *  $}• 

Note,  that  in  the  univariate  case  (15)  we  may  take  «  (tv,ty+^)  and  clearly  #J^ 

-  #Jv  -  k+1. 


Theorem  5.  Let  Z  be  a  regular , perturbation  of  a  triangulation  ( X,I ).  We  define  the  re¬ 
normalized  B-spline  basis  for  ^ f x, I 

m‘p,(*|Z)  -  M}MZ), 

(41) 

IM*P)I  -  1 
i  P 


and  assume  that  (c.f.  (39),  (40)) 


(42 ) 

Then  there  exists  a  constant  k  - 
-  s 

(43)  K  <  l  |a|p)1/p< 

s  j£lAk  3 


*  Jv  *  d  • 


c(m,k,p,K)  such  that 


»  l 

jslAk 


a  ,M*p) I  <  d1/p' 
j  )  P 


<  I  l-.lV/P. 

j£IAk  3 


Por  univariate  spline  functions  this  result  is  due  to  C.  de  Boor  [1-3]  and  has  become 
a  basic  tool  in  spline  approximation.  The  following  example  shows  that  the  assumptions  of 
Theorems  4,5  do  not  necessarily  imply  that  the  B-splines  M j ,  j  e  IAk,  correspond  to  a 
triangulation  of  Rm  x  A^  . 
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Example  6.  Consider  the  following  set  of  knots  which  can  be  extended  to  a  regular 
perturbation  of  a  triangulation  of  R2. 


x0  “  X1  ”  “  (-1,1), 

«  x^  -  (1,1),  x|  -  (l,e),  x*  «  (— 1 ,— e) 

The  balls  can  be  chosen  in  the  sets  ,  02  and  the  coefficient  of  the  B-spline 

M( *lxg«x^,xj ,x^)  in  the  expansion 


(16*) 

is  given  by  (c.f.  Theorem  3) 


'  JC)"! 


0 

0 

1 

1 

0 

0 

1 

1 

-1 

1 

1 

-1 

-1 

1 

1 

-1 

sgn  det 

-1 

-1 

-1 

1 

det 

-1 

-1 

e 

-E 

1 

1 

1 

1 

1 

1 

1 

1 

Hence  the  coefficient  may  vanish  or  become  negative  although  the  assumptions  in  Theorems 

4,5  are  satisfied.  Since  the  B-splines  are  linearly  independent  the  coefficients  in  (16’) 

2 

are  unique  and  therefore  there  cannot  exist  an  associated  triangulation  of  R  x  A  for 


e  >  0.  This  is  also  the  reason  why  we  cannot  use  the  coefficients  to  normalize  the  b- 


splines,  l.e. 


“  (0^ )  ,  as  in  the  univariate  case,  unless  we  make  some  more 


restrictive  assumptions  on  the  distribulon  of  the  knots. 


Proof  of  Theorem  5. 


By  (42)  we  have  (c.f.  (39),  (40)) 


,r  .  „<P>,P  <  ap/P'  V  i.  |P 

'1  *3  ”j  "o,U  4  d  l  '“j1 

V 


«p  „ 

3  p.u 
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Summing  this  inequality  with  respect  to  v  and  interchanging  the  order  of  summation  yields 


,  £a  M<PV<dp/p'  I  I  la  |p.«;p',p  -dp'p-  I  |a  |p. 

3  3  p  jdAk  v£jj  3  3  p  v  jelAk  3 

We  use  the  functionals  constructed  above  (c.f.  (36),  (38))  and  take  into  account  the 


|P.u'P> „P 


I  la  j  I 


different  normalization  of  the  B-splines 

,(p)  ,.(p),-1. 

Aj  *-  (  y  1  Aj  • 

(  D  ) 

We  shall  show  that  IA^  1^,  is  uniformly  bounded  which  implies  the  local  estimates 

>V  <  c  'S*3  *\P)'p,h  •  j£iik  • 

Summing  this  inequality  with  respect  to  j  e  IAk  gives  the  lower  estimate. 

Since  It  I  ,  (  c  we  have  by  (38) 
a  p’ 


.A<P,I  ,  <  c  (S'?1)-1  l  h7"/P-l“l 

j  P  j  l.foc1 


D  C  (OlX-x^  . 


We  first  observe  that 


(s^f1  <  c  h;m/p'. 


To  see  this,  consider  the  set  of  knots  {x  n}  of  the  B-spline  M . •  This  set  must 

r0  rn  3 

contain  points  x0,***,?"  such  that  0^  n  a(x°,**",xm)  *  t  <  since  otherwise 

fl  n  supp  M  «  *S  contradicting  (25).  (26)  implies  S  .  =  P  .  in  particular 

1  j  K I  K  | 

M  .  €  p  ,  Hence  we  must  have  J)  c  o(x° ,  •  •  •  jX*")  since  M.  is  not  smooth  across  the 

J  I  X  |  1  J 

(m-1)  dimensional  faces  of  this  simplex.  We  associate  with  11^  the  sinplex 

~  , ,  ~0  .  ,  -m  .  »1  ,«k 

o((x  ,eQ), •••,(><  ,eQ),(x  ,e1>,«»*,(x  ,e^)) 


where  {x  , ***,x  }  *  {x  ,•••,*  }  \  {x  , •••,x  }.  Clearly 

0  ~on 

vol  a  -  vol  o((x  ,eQ),***,  (xm,eQ ),( 0 ,  e1 ),•••,  (O.e^) )  -  h™  . 


and  by  Definition  1  we  obtain 


IM.I  <ch;m  h“/P  =  c  h-m/p' 
IP  i  i  i 


which  implies  the  estimate  (45). 


Now  consider  the  term  D  C,(0|X-x,).  By  its  definition  (c.f.  Theorem  3)  and  (22)  it 

3  i 


is  the  sum  of  terms  of  the  form 


m 


k 


where  Ib^l,  |z  I  <  <  .  This  yields  the  estimate 


(46) 


|DaCi(0|X-xi)l  <  c(ra,lt)(ich1)®f|al 


"j 


Combining  the  estimates  (45),  (46)  we  obtain 
(47) 


IA*P*  lpl  <  c(ra,)c,p)iem+l!  . 


Corollary  3.  Let  x,I  be  *  apace  of  splines  not  necessarily  satisfying  conditions 


(25)-(27).  Assume  that  for  a  family  of  bounded  functionals  L.:L  *  R,  supp  L  c  supp  M  , 

3  P  3  3 


IL . I  <  c  ,  the  associated  "quasiinterpolant” 

3  p 


(48) 


Qf  *-  l  (L  f'M'P) 
jclik  J  1 
reproduces  polynomials  of  total  degree  4  r,  i.e. 


(49) 

Then  we  have 

(50) 


Qp  -  p  ,  p  £  Pr. 


,Qf  -  f,p  fl  <  (dc+1)  di.tLp{5)(f,Pr) 


where  0  :«  u{supp  M^lsupp  n  0  *  $}  and  d  is  an  upper  bound  for  (c.f.  (39), 


(40)). 

Proof .  Denote  by  p  the  best  approximation  to  f  on  [1  .  By  (49)  we  have 


IQf-fl  <  IQ(f-p)l  +  lf-pl  . 
P  P  P 


The  first  term  can  be  estimated  by  (c.f.  (39),  (40)) 


■Q< 


"p"p  ■  l  '.l  4 1  a",r\L 1  J,  °p,,-p,Ulm)p,'p.v 


V  j«J 


v  V  jSJv  p€J^ 


Interchanging  the  order  of  summation  we  obtain  by  (41) 


i  j  i  i 

j  U€J  V«J  P'  It  J  P’  V  J  u£J  P'  U 

U  CI-+J.  J  J  J 


supp 


I  I  I  c'.f-P.J,,  . 

w  3£j„  u  u  u 
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V 


Corollary  3  is  the  essential  error  estimate  for  local  approximation  schemes  (c.f. 
[1/4]  for  the  univariate  and  [18]  for  the  tensor  product  case).  The  dual  functionals 
are  one  particular  choice  for  Lj  .  But,  they  are  merely  of  theoretical 
Importance.  Efficient  schemes  should  be  based  on  poin*  evaluation,  etc.  Of  course,  in 
principal,  such  methods  could  be  constructed  along  the  same  lines.  In  order  to  get  a 
smooth  approximation  the  major  problem,  however,  is  to  choose  the  knots  so  that  the 
derivatives  of  the  B-splines  are  not  too  big.  Moreover  there  should  be  efficient  ways  for 
computing  the  values  of  Ljf,  Mj  and  precise  estimates  for  the  constants  involved.  These 
problems  however  exceed  the  scope  of  this  paper  and  we  give  only  one  example  illustrating 
the  above  concept.  The  scheme  constructed  below  is  in  the  univariate  case  due  to  M.  J. 
Marsden  and  I.  J.  Schoenberg  [20]  and  for  two  variables  to  T.  N.  T.  Goodman  and  S*  L.  Lee 
[15). 


Example  7.  Consider  the  following  special  cases  of  identity  ( 1 9 1 >  (c.f.  also  Theorem  3  and 
(21)) 

1  5  I  C.  M 


jelAk 


J  j 


(51) 


and  let 


'xL""’  v*1 


£j  -  cj’  k'’ 


We  define  the  approximation  scheme 

(52) 

which  by  (51)  reproduces  linear  functions.  If  we  furthermore  assume 

(53)  0^  >  0  ,  5^  £  supp 

then  Q  is  a  positive  linear  operator  and  we  obtain  by  a  slight  modification  of  Corollary 
3 


Of  -  I  fly  cj  Mj 


(54) 


where  h  sup  dlam  supp  M 

j 


IQf-fl  <  c(m,k)  h  |f| 


2,« 


J 
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In  the  univariate  case,  condltlona  (S3)  ara  automatically  satisfied.  in  general  we 


obtain  from  (51),  if  the  B-splines  aca  linearly  independent, 

VT  “  cj_1  k’1  ■fcucj<y,x-t,lyo  • 

Let  t  be  the  canter  of  supp  Mj  «  By  (46)  it  la  therefore  neceaaary  for  e  eupp 
that 


<  c  (diam  supp  M^) 

which  leads  to  restrictions  for  the  distribution  of  the  knots  (c.f,  Exaaple  6). 
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